Effective models for gapped phases of strongly correlated quan- 
tum lattice models 
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Abstract. - We present a robust scheme to derive effective models non-perturbatively for quan- 
tum lattice models when at least one degree of freedom is gapped. A combination of graph theory 
and the method of continuous unitary transformations (gCUTs) is shown to efficiently capture 
all zero-temperature fluctuations in a controlled spatial range. The gCUT can be used either 
for effective quasi-particle descriptions or for effective low-energy descriptions in case of infinitely 
degenerate subspaces. We illustrate the method for Id and 2d lattice models yielding convinc- 
ing results in the thermodynamic limit. We find that the recently discovered spin liquid in the 
Hubbard model on the honeycomb lattice lies outside the perturbative strong-coupling regime. 
Various extensions and perspectives of the gCUT are discussed. 



Introduction. — A microscopic description of 
strongly correlated quantum many-body systems is very 
demanding and powerful theoretical techniques are clearly 
needed. Especially in two dimensions there is a lack of the- 
oretical tools to tackle microscopic models. Despite quan- 
tum Monte Carlo, functional renormalization group, and 
2d extensions of density matrix renormalization group, 
the method of continuous unitary transformations (CUTs) 
[1-4] has been applied successfully to a variety of physi- 
cal systems. Here a CUT is used to derive effective low- 
energy models with the major task to solve an infinite 
hierarchy of differential equations. One controlled way to 
proceed is to truncate these equations in a perturbative 
manner yielding a high-order series expansion of the effec- 
tive low-energy model [5-7]. Such a perturbative version 
(pCUTs) is similar in spirit to linked cluster expansions 
(LCEs) which have become a standard tool to treat a large 
class of quantum lattice models [8,9]. LCEs exploit the 
linked-cluster theorem by calculating on finite graphs the 
correct amplitudes of the effective low-energy model in 
the thermodynamic limit up to high orders in perturba- 
tion. In contrast to conventional LCEs, the pCUT has 
been also applied succesfully to infinitely degenerated sub- 
spaces [10-12] and to topologically ordered ground states 
[13]. 

The applicability of pCUTs is clearly limited by its per- 
turbative nature. A non-perturbativc extension is there- 



fore highly desirable. Attemps in the past solved the flow 
equations in a self-similar way (sCUTs) by truncating in 
terms of the spatial extension and the structure of the 
operators [14-17]. The major complication with these ap- 
proaches is to control the error induced by truncations. 

Here we present a graph-theory based continuous uni- 
tary transformation (gCUT) which yields a robust scheme 
to derive effective models non-perturbatively by integrat- 
ing out gapped degrees of freedom in a controlled spatial 
range. It can be applied for all kinds of lattice topologies, 
for all kinds of interacting degrees of freedom, and for de- 
riving low-energy models with a reduced Hilbert space. 

Let H be an Hamiltonian on a given arbitrary lattice at 
zero temperature. On each vertex of the lattice a quantum 
degree of freedom like a spin, a boson, or a fermion is 
placed which interacts via short-range operators. 

We use a CUT to map H unitarily to an effective Hamil- 
tonian H c g. In most cases we focus on a quasi-particle 
(QP) conserving CUT [5,6], i.e. the effective Hamiltonian 
is blockdiagonal in a QP counting operator Q. Formally, 
we can write H e s = J2^Lo -^eff > wnere Q denotes the num- 
ber of QPs. The low-energy physics is then expected to 
be contained in the blocks with q small simplifying the 
many-body problem tremendously 

Each QP block of the effective Hamiltonian H^g = 
J2i a ; 9 c ff^'i 9 ' ) 1S generically a sum of g-particic operators 
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weighted by an amplitude afl s . The final goal of 
the gCUT is to calculate these amplitudes quantitatively 
in the thermodynamic limit. The gCUT is defined by the 
following steps: 

(i) Generate all topologically distinct connected graphs 
Q v on the lattice and sort them by their number of sites 
n. 

(ii) For each graph Q v one sets up the finite number of 
flow equations using the method of CUTs [1-4]. A con- 
tinuous auxiliary variable 1 is introduced defining the 1- 
dependent Hamiltonian H G » (l):=W (l)H g »U(l). Then the 
flow equation is given by diH g "(l) = [fj(l), H g " (/)], where 
fj(l):=—U'(l)(diU(l)) is the anti-Hermitian generator of 
the unitary transformation. At the end of the flow one ob- 
tains an effective graph-dependent Hamiltonian ff% with 
amplitudes a^g. For the QP conserving generator, H^X 
commutes with a chosen QP counting operator Q and the 
amplitudes aSp can be sorted by the number of QPs in- 
volved. In an eigenbasis of the operator Q, the QP con- 
serving generator is given by fj^(l) — sgn (g, — qj) Hij(l), 

where qi and qj are eigenvalues of Q. One can also 
use the generator introduced by Wegner [1] which reads 
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where the diagonal part 



H d is given by all matrix elements between states with 
the same number of QPs and H nd denote all remaining 
non-diagonal processes. Since no truncations are done for 
a given graph, the flow equations converge [1,14,18] and, 
consequently, the gCUT is intrinsically robust. 

(iii) The graph-dependant amplitudes contain con- 
tributions from smaller subgraphs with fi < v which 
have to be substracted in order to avoid double counting. 
We therefore define reduced amplitudes a^{v) represent- 
ing all fluctuations where the whole graph with all n sites 
are involved 



(1) 



The next step is to collect properly all the reduced con- 
tributions from the different graphs in order to get the 
amplitudes a!f2 in the thermodynamic limit. One finds 



ag = £ o$ (») = £ E m(u) og (u) • (2) 

n=l n=l \u\=n 

Here m{v) denotes the integer number how many times 
the process belonging to the amplitude (v) on graph v 
can be embedded on the real lattice. 

The demand of the full gCUT calculation scales mainly 
with the maximal number of sites n max of the largest 
graphs and their total number. The number of flow equa- 
tions to solve is given by the square of the symmetry- 
reduced Hilbert space on a given graph. The number n max 
corresponds physically to the spatial range of fluctuations 
treated in the gCUT calculation. So as long as the corre- 
lation length is finite which is true for any gapped system, 



the gCUT enables a quantitative and non-perturbative 
derivation of effective low-energy models when n max is cho- 
sen large enough. 

Before discussing explicitly applications of the gCUT, 
we present some generic and important properties. In gen- 
eral, one should distinguish two types of applications: i) 
QPs above a unique ground state and ii) effective low- 
energy models for infinitely degenerated subspaces. 

In the first case, one is interested in the QP proper- 
ties above a strongly correlated ground state. Here the 
ground state is up to topological degeneracies unique and 
it can be characterized by the ground state energy per 
site which is given by the q = block of H c s- The one- 
particle dispersion can be obtained from the q = 1 block 
and, consequently, multi-particle properties are contained 
in blocks with q > 1. 

Let us mention that for the single-valued q = block 
the gCUT just diagonalizes the problem and the calcula- 
tion then corresponds to exact diagonalization of the finite 
graphs. Such an exact linked cluster expansion for the 
ground state energy has been already discovered and used 
successfully by Irving and Hamer in 1984 in the context 
of lattice gauge theories [19], and it has been extended 
to thermodynamic quantities for quantum lattice models 
recently [20]. But the gCUT goes well beyond these tech- 
niques for q > 0. 

In the second case, the system consists of different de- 
grees of freedom and only some of them are gapped. Then 
one is interested in deriving an effective model in the re- 
duced but still infinitely large q = Hilbert space. One 
typical example which we will detail below is the strong 
coupling limit of single-band Hubbard models where an ef- 
fective description in terms of interacting quantum spins 
is appropriate since charge degrees of freedom are gapped. 

It is important to understand the difference of the gen- 
erators ?y qp and 7y We g ner . Let us distinguish for both gener- 
ators the states of the full Hilbert space on a certain graph 
according to the unperturbed problem: We define the low- 
est energy states to have q — and all the remaining states 
to have q > 0. One then has to distinguish two regimes: a) 
If all energies corresponding to the unitarily transformed 
q = states are energetically below the rest of the states, 
then both generators give successfully an effective low- 
energy description. The effective models obtained by the 
two generators can differ by a unitary transformation of 
the q = Hilbert space but yield the exact spectrum of 
the considered graph, b) In contrast, if the spectrum of a 
graph is such that q = states and q > states overlap 
energetically, then the two generators behave differently. 
It is well known for r? qp [5] that it sorts the eigenvalues 
in ascending order. Then q = states with a high-energy 
are replaced by low-energy q > states in the low-energy 
block during the solution of the flow equation. As a con- 
sequence, no well-defined effective low-energy theory can 
be derived with ?y qp . This is different for 7y w °g ncr which al- 
ways keeps track of the unitarily transformed q = states 
and one is able to deduce an effective low-energy model. 
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The different behaviour of both generators is explicitely 
found for the Hubbard model on the honeycomb lattice 
(see below). 

The gCUT for q = is similar to the contractor renor- 
malization group (CORE) [21]. Indeed, in CORE the ex- 
act low-energy spectra of graphs for the full model are 
matched to the low-energy spectra of an effective Hamil- 
tonian in a reduced Hilbert space. The obtained effec- 
tive model is identical to the one obtained by gCUT up 
to a unitary transformation inside the low-energy Hilbert 
space. Interestingly, also in CORE one has to be very 
careful if states classified as low- and high-energy states in 
a perturbative limit overlap energetically [21]. If one just 
chooses the states having the lowest energy on a given 
cluster, then no well-defined effective model can be found 
like with 77 qp . So also in CORE one has to find a way 
to choose the relevant states in order to build an effective 
low-energy model [21]. This is actually done automatically 
in the gCUT approach when using 7y We s ner . 

The gCUT can be implemented for any generator fj 
which is optimal for a given physical problem. If one is 
for example only interested in certain QP blocks or would 
like to treat the decay of QPs, it might be advantageous 
to only partially block diagonalizing the problem. This 
is exactly achieved with the generators proposed recently 
[16,17]. 

As a first example for type i), we apply the gCUT for 
the Id Ising model in a transverse field which is given by 



H 



TFIM 



(3) 



with J assumed to be positive. This model is exactly solv- 
able by fermionization [22] and it allows therefore a quan- 
titative gauge of the gCUT. Additionally, it nicely serves 
as a pedagogical example where the elementary proper- 
ties of the gCUT can be demonstrated explicitely. The 
model is self-dual and one has a quantum phase transi- 
tion at J = 1. It is therefore sufficient to concentrate on 
the parameter space J < 1. The unique ground state is 
unitarily connected to the polarized state where all spins 
point in z-direction and elementary excitations are dressed 
spin flips. Consequently, the operator Q is chosen to be a 
counting operator of down spins. 

The graph theory aspects are particularly simple in Id. 
Graphs arc chain segments and there exist only one graph 
for a given n (see Fig. la). All graphs v can be embedded 
JV-times on the lattice and one has miv) — N. We have 
calculated the OQP and 1QP sector up to n max = 11. 

The q = sector is a single number which corresponds 
to the ground state energy per site e/N. Following the 
gCUT recipe, one has for a chosen n max 
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Fig. 1: (Color online) (a) Illustration of the first three graphs 
with n G {1, 2, 3} for the Id Ising model in a transverse field. 
The sites on the graphs are denoted by a, /3, 7, et cetera from 
left to right, (b) The local hopping amplitude ao as a functions 
of the graph size n is displayed for J = 0.9. The gCUT results 
are denoted by triangles and the exact value is indicated by 
the solid line. The inset shows the difference between the exact 
value and the gCUT results on a logarithmic scale. 

where Effi is the ground state energy of graph v obtained 
by the gCUT and e^(u) is the corrected energy where 
contributions from all subgraphs are substracted. For 
n max < 3 one finds the following expressions 
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The corresponding results up to n max = 11 arc shown in 
the upper panel of Fig. 2. The gCUT results monotoni- 
cally improve for increasing n max . As expected, the largest 
deviations are found at the critical point J = 1. But for 
J < 1, the gCUT converges quickly to the exact result. 

The q = 1 sector contains the one-particle dynam- 
ics in the form of hopping elements ag. The ampli- 
tude as corresponds to a hopping from site i to site 
i + S on the infinite lattice. A gCUT calculation for a 
given n max includes therefore hopping elements as with 
S € {0, ±1, . . . , ±(n max — 1)} and one finds 
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Fig. 2: (Color online) Comparison of exact and gCUT results 
for the Id Ising model in a transverse field: (a) ground state 
energy per site as a function of J. The inset displays differences 
between exact and gCUT results for various n max . (b) One- 
particle dispersion cj(k) for different values of J. Solid lines 
represent the exact results while symbols correspond to gCUT 
results. Upper inset represents a zoom. Lower inset shows the 
one-particle gap A = u(0) at J = 0.9. 

where a^j +& e s{ v ) denotes the effective hopping amplitude 
from site i to site i + S on the graph v where all subgraph 
contributions are already substracted. 

We would like to illustrate this procedure for the local 
hopping ciq and for n max < 3. We denote the sites on the 
graphs by a, f3, 7 et cetera from left to right (sec Fig. 2a). 
The case n max = 1 is then simply given by 

a (1) Jl) = a Gl „ . 

a,CK,effV / a, a, erf 

For n max = 2 one finds for the two amplitudes the follow- 
ing expressions 

a (1) J2) = a G2 „-a {1) Jl) 

a, a, eft V / a, a, eft a,a,crtv / 

a ie,^,eff (2) = a %fi,aS ~ a a ,a,eS ("*-) ' 

In practise, we stress that there are of course symmetries of 
the problem on a given graph which should be considered 



in order to reduce the effort of the calculation. In the 
present case all chain segments have a reflection symmetry 
with respect to the geometrical center of the chain. This 
immediately implies a^ aeS (2) = a^ eff (2). Finally, the 
analogue formulas for the three processes on the graph 
with three sites read 

a i,Leff (3) = a Q^a,eff — a La,cff (2) ~ a i,Lcff W 

(1) /q\ Gz 

a /3,/3,effW — a /3,^,eff 

a (1) M) = a Gi „-a [ n\ J2) - a (1) Jl) . 

The corresponding results for the local hopping element 
ao for graphs up to n max = 11 are displayed in Fig. 2b 
and compared to the exactly known value. 

The one-particle dispersion uj(k) is obtained by Fourier 
transformation u(k) = ao + 2 X^fci a <5 cos(<5fc). The gCUT 
results are displayed in the lower panel of Fig. 2. The 
gCUT can again be compared with the analytically known 
exact result 2i/l + J 2 — 2Jcos(fc). The only sizable de- 
viations occur close to the critical point J = 1 and close 
to k = 0. Otherwise, the gCUT one-particle dispersion 
is basically converged and captures all relevant aspects of 
the exact expression. 

As an example of type ii) , we study the Mott phase of 
the single-band Hubbard model on the triangular and the 
honeycomb lattice at half filling. The model is given by 

& = u n 't"4 - * ( c t c J° + h - c -) • ( 6 ) 

In the strong coupling limit t/U — > 0, the system has a 
finite charge gap and it is expected that the physics can 
be described solely by spin degrees of freedom [23,24]. 
Indeed, leading second order perturbation theory gives a 
nearest-neighbor Heisenberg model. But it is an inter- 
esting and a challenging problem whether one can derive 
an effective spin model also close to the metal-insulator 
transition when the Mott insulator is already weak and 
corrections to the Heisenberg exchange become important. 
Physically, interesting spin liquid (SL) phases are expected 
to be realized in this regime [12,25-27]. 

The operator Q is chosen as a counting operator of 
double occupancies. In contrast to the previous example 
where the main interest was in blocks with a finite num- 
ber of QPs, here we focus on the block with no double 
occupancies which can be expressed as a spin 1/2 model 

-^eff = J2i,H Jn (Si ■ S i +f ^j +. . ., where i and i+n a denote 
sites on the lattice. All remaining terms can be written 
in an analog way as products of ■ Sj due to the SU(2) 
symmetry of the Hubbard model [12, 15]. 

In the following we focus on the nearest neighbor and 
next-nearest neighbor Heisenberg exchanges J\ and J-2- 
We calculated both couplings for the two lattices with 
gCUTs up to n max = 6 using the qp and the Wegner 
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generator. As explained above, the generators behave dif- 
ferently if unitarily transformed q = and q > states 
overlap energetically. We have found this to be the case 
for gCUT(6) around t/U = 0.25. In this range of param- 
eters one has to use the Wegner generator to obtain an 
effective spin model. In all other cases we studied the re- 
sults obtained by both generators are very similar. In total 
there are 67 (11) different graphs for the triangular (hon- 
eycomb) lattice. Additionally, we have performed a pCUT 
calculation up to order 12 (14) in t/U for the triangular 
(honeycomb) lattice. The results are shown in Figs. 3-4. 

We remark that the gCUT curves for different n max 
smoothly converge in a wide range of t/U parameters. In- 
terestingly, the bare pCUT series become unreliable well 
before and it is mandatory to use extrapolation schemes. 
We have found that different approximations like Pade or 
DlogPade work well for the nearest neighbor exchange J\ . 
In contrast, such standard extrapolation schemes fail for 
the next- nearest neighbor exchange Ji in the t/U range 
of the SL on the honeycomb lattice due to the appear- 
ance of spurious poles. We therefore have applied so-called 
self-similar appromixants [28] . A self-similar approximant 
SS(fc) of order k is defined as the function 



P 0.05 ■ 



SS(k) = (1 + Aix) ni (1 + A 2 x) n * . . . (1 + A k ) r 



(7) 



where Ai and rz.j are complex coefficients. These coeffi- 
cients have to be determined by comparing the Taylor ex- 
pansion of SS(k) in the variable x to the obtained pCUT 
series with x = t/U. We found convincing and convergent 
results for both couplings on both lattices. The largest de- 
viation between two different subsequent self-similar ap- 
proximants SS(k) and SS(k + 1) are found for J 2 on the 
honeycomb lattice as displayed in Fig. 4b. Remarkably, 
the gCUT and the self-similar approximants are in very 
good agreement for both couplings on both lattices for 
the range of considered t/U values. 

This is a very promising result, since it shows that the 
gCUT can be at least as good as extrapolations of high- 
order pCUT series. For the honeycomb lattice, we stress 
that the recently discovered SL [27] lies not in the purely 
perturbative regime in contrast to the one on the trian- 
gular lattice [12]. One attractive scenario put forward in 
the literature [29,30] is the relevance of a 3\-Ji model for 
the SL phase on the honeycomb lattice. It is argued that 
due to the increasing frustration when increasing t/U the 
long-range ordered Neel state found at small values oit/U 
is unstable towards a SL phase. We find J2/J1 ~ 0.06 at 
the critical point U/t — 4.3 obtained from quantum Monte 
Carlo simulations [27]. It should be clarified in the future 
whether this rather small ratio is sufficient to destabilize 
the Neel order on the honeycomb lattice or whether other 
spin-couplings are relevant for the realization of the SL. 

Both examples show that the gCUT is a very robust 
and flexible tool to derive effective low-energy models in 
a quantitative fashion. The gCUT is solely limited by the 
maximal n max which can still be treated. One therefore 
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Fig. 3: (Color online) Nearest-neighbor Heisenberg exchange J 1 
as a function of t/U for the (a) triangular lattice and (b) the 
honeycomb lattice. We always compare bare (black lines) and 
extrapolated order n self-similar approximants of the pCUT 
series SS(n) to gCUT curves with varying n max using the qp 
(filled symbols) and the Wegner generator (open symbols). The 
first (second) dashed vertical line illustrates for both lattices 
the beginning (end) of the SL [12,27]. 



demands that the correlation length is not larger than the 
range of fluctuations captured by the treated graphs. 

If larger n max are needed to capture the physics, one has 
to develop strategies to treat larger n. In practise, this is 
most likely possible by reducing the number of flow equa- 
tions such that on the one hand the graph becomes nu- 
merically tractable and on the other hand the amplitudes 
of interest are still quantitatively correct. One therefore 
has to introduce truncations which are most physically 
performed like in current sCUT calculations where blocks 
belonging to a large number of QPs are discarded from 
the calculation because their influence on the blocks with 
a few number of QPs is expected to be small. But we note 
that the current gCUT set up has some intrinsic advan- 
tages: i) The existence of a fix point for the full graph is 
guaranteed by the finite size. One can therefore stop the 
flow once a convergence of the diagonal couplings of inter- 
est is achieved, ii) The quality of these amplitudes can be 
checked quantitatively by using exact diagonalization for 
the full model (as long as n is not too large). 

It is one strength of CUTs that also effective observ- 
ables and dynamical correlation functions can be calcu- 
lated. To this end any observable has to be transformed 
with the same unitary transformation as the Hamiltonian. 
A gCUT for observables amounts therefore to perform the 
same steps as for the Hamiltonian: i) collect all graphs 
with fixed n. ii) Solve flow equations for the observable 
on the graphs, iii) Do subcluster substractions and em- 
beddings on the real lattice. 

There are actually many physical problems where gCUT 
is expected to give important new insights: 
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Fig. 4: (Color online) Next-nearest neighbor Heisenberg ex- 
change J2 as a function of t/U for the (a) triangular lattice 
and (b) the honeycomb lattice. Notations are equivalent to 
the ones used in Fig. 3. 



1) The pCUT needs an equidistant unperturbed spec- 
trum. The gCUT does not have this limitation. A much 
larger class of microscopic problems can therefore be tack- 
led including quantum spin models with S>l/2, models of 
mixed degrees of freedom like the t-J model giving rise to 
unconventional superconductivity, multi-orbital Hubbard 
models with more than one local energy scale U important 
for realistic modelling of materials, or bosonic models rel- 
evant for cold atoms in optical lattices. 

2) The gCUT allows also the calculation of dynamical 
correlation functions for a much larger class of problems. 
This opens the possibility to obtain a microscopic under- 
standing of various spectroscopic experiments like it has 
been done successfully for many pCUT applications. 

3) The gCUT can also be applied in situation where no 
obvious perturbative limit (as for the discussed examples) 
can be found. It is only important that an appropriate 
counting operator Q and its eigenbasis can be defined. 
Note that also symmetry broken phases can be tackled. 

4) The gCUT is especially well suited also in three di- 
mensions since n max and therefore the radius of fluctua- 
tions kept in the gCUT is basically constant since only the 
number of graphs for a given n max increase. 

In conclusion, we developped a robust scheme to de- 
rive non-perturbatively effective Hamiltonians for strongly 
correlated quantum lattice models. Such gCUTs combine 
graph theory and CUTs, and can be applied to any lattice, 
in any dimension, and for any degree of freedom as long 
as one sort is gapped. It can be either used to study the 
properties of elementary excitations of a strongly corre- 
lated system or to derive effective low-energy Hamiltoni- 
ans for infinitely degenerated subspaces. The latter might 
then be tackled with other powerful techniques like cluster 
dynamical mean-field theory [31] or functional renormal- 
ization group [32,33]. 
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